Shape changes and large-amphtude collective dynamics in neutron-rich Cr isotopes 
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Shape-phase transition in neutron-rich Cr isotopes around A'^ = 40 is investigated by employing the 
collective Hamiltonian approach. The inertial functions for large-amplitude vibration and rotation 
are evaluated by the local normal modes along the axial quadrupole collective coordinate using the 
Skyrme and pairing energy density functionals. The time-odd components of the mean fields are 
fully included in the derived masses. The low- lying spectra obtained by requantizing the collective 
Hamiltonian show an excellent agreement with the recent experimental data. Gradual change from 
spherical to axially deformed shapes in between A'' = 34 and 38 is well described. 
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PACS numbers: 21.10.Rc; 21.60.Ev; 21.60.Jz; 21.10.Ky 

How can we define the shape of quantum many-body 
systems? Atomic nuclei have a variety of equilibrium 
shapes, and show shape changes along the isotopic or 
isotonic chain. The mean-field approximation gives us 
an intuitive picture of the nuclear shape. However, we 
have to go beyond it to describe the shape-phase tran- 
sition; the dynamical change of the mean-field potential 
associated with the large-amplitude collective motion. 

The self-consistent mean-field model employing the 
effective interaction or the nuclear energy-density- 
functional (EDF) method has successfully described the 
ground-state properties Recent advances in the 

computing capability together with the highly-developed 
techniques in the nuclear EDF method allow us to calcu- 
late the ground-state properties of nuclei including defor- 
mation in the entire mass region of the nuclear chart Q ■ 

Magic number or shell closure is an essential concept 
in understanding the stability against the deformation. 
The subshell closure at 40 created by the gap between 
1(79/2 and 2pi/2i I/5/2 orbitals has attracted much atten- 
tion for several reasons 0]. The proton-rich N = Z nu- 
cleus ^°Zr lies in the center of the well-deformed ^ ~ 80 
region . This is because a shell gap of 40 again appears 
in a deformed region, and the shell effect of protons and 
neutrons coherently stabilize the nucleus deformed. On 
the other hand, the existence of -/V = 40 subshell closure 
is suggested experimentally for the neutron-rich nucleus 
^^Ni [sl, 0] . The strength of this subshell closure and its 
persistence for Z < 28 determine the waiting point for 
the r-process nucleosynthesis at ^^Cr, which is considered 
to be a progenitor of ^^Ni 0] . 

The half-life measurement at CERN-ISOLDE have de- 
duced that the neutron-rich ®^Fe is deformed with a 
quadrupole deformation /? ~ 0.26 Q. Since the Cr iso- 
topes lie at mid proton I/7/2 shell, protons could addi- 
tionally destabilize the nucleus in a spherical shape and 
favor deformation. Experimental evidences of the nuclear 
shape changes are related to low-lying quadrupole collec- 
tivity, such as the small excitation energy of the 2^ state, 
the ratio of the excitation energy of the 2^ and 4^ states 



B[E2;2^ 0^), etc. The observed small excitation en- 
ergy of the 2^ state in neutron-rich Cr isotopes indicates 
that the deformation develops toward = 40 

Very recently, the large-scale shell-model calculation 
for the neutron-rich Cr isotopes has become avail- 
able |14| . where the Hamiltonian with dimension of ^ 
10^° is diagonalizcd. However, deducing the picture of 
deformation is difficult in the fully quantum-mechanical 
calculations. 

In this article, on the basis of the nuclear EDF 
method closely related to the mean-field approxima- 
tion, we develop a new framework of the microscopic 
theory for the large-amplitude collective motion em- 
ploying the quadrupole collective Hamiltonian approach 
for the axially-symmetric nuclei. The similar attempts 
of the EDF-based collective Hamiltonian starting from 
the Gogny interaction 15, [l^, the relativistic La- 
and the Skyrme interaction [l9| have 
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been made for description of the large-amplitude col- 
lective motion. The cranking approximation, however, 
has been applied to calculate the inertial functions 
[lij . and the time-odd components of the mean field re- 
main largely unexplored except for a few attempts in the 
adiabatic time-dependent Hartree-Fock-Bogoliubov the- 
ories Our method includes the time-odd mean 
fields in the inertial functions. And this method is ap- 
plied to the shape-phase transition in neutron-rich Cr 
isotopes around A^ = 40. 

We start from the (1 + 2)D (vibration along the (3 di- 
rection and rotations about the two axes which are per- 
pendicular to the symmetry axis) quadrupole collective 
Hamiltonian; 
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i^J,(/3Vf-fn/3)- (1) 
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-E4-1-/-E2+, and the reduced transition probability 



With the Pauli's prescription, we requantize Eq. ([T]) and 
construct the collective Schrodinger equation; 

HcoM^ c.Im{P, 91,02) ^ Ec.I-^a:IM{P, 91,62), (2) 

where ^aiM{P, ^1,^2) is the collective wave functions in 
the laboratory frame, / is the angular momentum quan- 
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turn number, M is the z— component of 7, and E^i is the 
excitation energy. The cohective wave functions are writ- 
ten in terms of wave functions in the body-fixed frame 

«'a/if(/3); 

aiK{P){^\IMK) = $Q/(/3)(r2|/MO), 

(3) 

with {Q,\IMK) a superposition of the rotational wave 
functions. 

The collective Schrodinger equation in the intrinsic 
frame now reads 

Tvib + y^li^ + V{P)\ = Ea,i<^^i{P), (4) 



2J(/3) 



where 
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2A/^(/3) a/32 



2Af^(/3) 



d_ 

(5) 



and the vibrational wave function is normalized as 



with the metric |G(/3)| = Mp{l3)J^{l3). 

The collective potential V{I3) is calculated by solving 
the constrained Hartrec-Fock-Bogoliubov (CHFB) equa- 
tion; 

s{m\HcKFB\m) = o, (7) 

-ffcHFB =-ff-5I^"^"-MQ20. (8) 



The microscopic Hamiltonian H is constructed from 
the Skyrme and pairing EDFs. Following the discussion 
in Ref. [23|, the vibrational mass M^(/3) is evaluated us- 
ing the local normal mode as 



(9) 



T 9Q20 0Q20 

^ = ■^m)\Q2,\m) = {m\[Q2o.^]\m) 

(10) 

by solving the local quasiparticlc-random-phasc approx- 
imation (LQRPA) equation on top of the CHFB state 

\m) ; 



\m) = o, (11) 



6{m\[HcHFB, —] - ujtQ,Am) = 0, (12) 



(13) 



Here, the isoscalar quadrupole moment operator Q20 s-nd 
1] = yJH I ^f^A^r^) are calculated at each state |0(/3)). We 
write the collective coordinate as q^. The rotational mo- 
ment of inertia »/(/?) is evaluated by the LQRPA equa- 
tion for collective rotation (extension of Thouless-Valatin 
equation to CHFB states [23[). 

The reduced quadrupole matrix elements used for cal- 
culating the E2 transition probability and the spectro- 
scopic quadrupole moment are evaluated as 



(a/||i?2||a'/') =v/(2/+l)(2/' + l)(-)^ 

7 2 
.0 



{(iI\Fe2\oII' 



(14) 



i,ocI\FE2\olI') = y d/3<i>;,(/3)^^B2(/3)$a'7'(/3)|G(/3)^/^ 

(15) 

where Fe^^P) = {4'{(^)\Pe2\4>{P)) j and Fe2 is the electric 
quadrupole moment operator. 

We solve the CHFB plus LQRPA equations employing 
the extended procedures of Rcfs. [25j, |2^]. To describe 
the nuclear deformation and the pairing correlations si- 
multaneously with a good account of the continuum, we 
solve the CHFB equations in the coordinate space using 
cylindrical coordinates r — (p, z,(j)) with a mesh size of 
Ap = Az = 0.5 fm and a box boundary condition at 
(Pmax, -Zmax) = (12.25, 12.00) fm. Wc assumc axial and 
reflection symmetries for the CHFB states. The quasi- 
particle (qp) states are truncated according to the qp 
energy cutoff at Ea < 60 MeV. We introduce the addi- 
tional truncation for the LQRPA calculation, in terms 
of the two-quasiparticle (2qp) energy as Ea + Ep < 60 
McV. In the present calculation, the LQRPA equations 
are solved on top of around 20 CHFB states. For each 
LQRPA calculation, we employ 256 CPUs and it takes 
about 185 CPU hours to calculate the inertial functions. 
Thus, for construction of the collective Hamiltonian, it 
takes about 3 700 CPU hours for each nucleus. 

For the normal (particle-hole) part of the EDF, we 
employ the SkM* functional , and for the pairing en- 
ergy, we adopt the volume-type pairing interaction as in 
Ref. [23|, where the deformation mechanism in neutron- 
rich Cr isotopes is investigated in detail. 

Figure [T] shows the total energy curves (collective po- 
tentials) of neutron-rich Cr isotopes under investigation. 
In ^*Cr and ^°Cr, the potentials are soft against the (3 
deformation even we get the HFB minima at /3 = 0. Be- 
yond = 38, the prolate minimum gradually develops 
up to A = 44. We can see a shoulder or a shallow mini- 
mum at the oblately deformed region as well. It is noted 
that the systematic calculation predict ^^Cr is spherical 
due to the magicity at A^ = 50 [3| using the same EDF 
as in the present calculation. 

Figure [2] shows the properties of low-lying states. 
When the neutron number increases from ^^Cr, the ex- 
citation energy of the 2f state drops toward ^^Cr. The 
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FIG. 1: Total energy curves of neutron-rich Cr isotopes cal- 
culated with SkM* and volume pairing interactions. 



i?4/2 ratio and the transition probability B{E2; 0^ — > 
2'^) grow up at the same time. Although we overestimate 
E{2^) at N = 34, these features are consistent with the 
experimental results [lll - [l3t , and it clearly shows the on- 
set of deformation at ~ 38. The experimental B{E2) 
value is available only in ^*Cr The present calcula- 
tion underestimates the observation. Underestimation of 
the quadrupole collectivity of ^^Cr probably attributes 
to the restriction of the collective coordinate to the /3 
and two rotational degree of freedom and / or to the EDF 
employed in the present calculation. 

Let us consider here the simplified case in solving 
Eq. (|4]) with a harmonic oscillator potential and a con- 
stant vibrational mass; 



B, J{/3) = 35/3^ 



(16) 



We obtain E{2+) = 2.00 McV and S(4+) = 4.26 McV, 
when we adopt B = 50 MeY'^ an d C = 200 MeV. The 
exact value for i;(2]^) is w = y^C/B = 2 MeV and 
£"(4^ ) — 2uj — A McV in the full fivc-dimensional col- 
lective Hamiltonian. This simple exercise implies that 
the present (1 -I- 2)D collective Hamiltonian can describe 
reasonably well the vibrational spectra of spherical nu- 
clei. However, the B{E2] 2f 0^) value in the (1 + 2)1) 
model is 3/5 times as large as an exact one. Therefore, in 
(nearly) spherical nuclei, all the five quadrupole degrees 
of freedom should be treated on the same footing, and it 
remains as future work. 
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FIG. 2: (Color online) (a): Excitation energies of the 2^ state 
in Cr isotopes, (b): Ratios of the £'(4^'') to E{2+). (c): Re- 
duced transition probabilities B{E2;0i — > 2^). (d): Spectro- 
scopic quadrupole moment of the 2f state, (e) ; Ratios of the 
transition matrix elements of neutrons to protons divided by 
the neutron to proton numbers, (f): Root-mean-square mat- 
ter radii. Experimental data are taken from Refs. [lOl - ilSll . 
Results of the shell-model calculations (denoted as SMI [2g] 
and SM2 [ij) are also included. 



The results of the spherical shell-model calculations 
employing the pairi ng-p lus- multipole forces with the 
monopole correction |28t and the LNPS interaction 
are also shown in Fig. [2l Lowering of E{2f) toward 
iV = 38 is consistent with the experiments and with the 
present calculation. 

We show in Fig. [2je) the ratios of the matrix elements 
of neutrons to protons M^/M^^ divided by the neutron to 
proton numbers for the excitation to the 2^ state. If the 
neutrons and protons coherently contribute to the exci- 
tation mode, M^/Mtj may approach N/Z. We can see 
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FIG. 3: (Color online) (a): Quadrupole masses calculated by 
the LQRPA and the cranking formula in ®*Cr. (b): Excitation 
energies squared for the associated collective modes. 



that the neutron excitation develops from = 34 to 38 
associated with the increase of the quadrupole collectiv- 
ity. 

Evolution of deformation is also seen in the matter 
radii as shown in Fig. [2][f). For the reference, we show 
the A^/^ dependence and the HFB results by a solid line 
and open circles, respectively. In the deformed systems, 
the increase of radius is proportional to Thus, beyond 
= 38 the radius is constantly larger than the systemat- 
ics of A^/"^. Furthermore, we can see that the calculated 
radius is larger than the HFB result in ^*Cr and ^"Cr 
associated with the large-amplitude fluctuation. 

It is interesting to compare our results with those 
in Ref. where the collective Hamiltonian approach 
is also employed. In Ref. they employ the five- 

dimensional collective Hamiltonian, in which both the 
P and 7 degrees of freedom are included. And the DIS- 
Gogny EDF is used for the particle-hole and particle- 
particle parts of the energy density. The rotational mass 
is evaluated by the Inglis-Belyaev procedure, and the vi- 
brational mass is calculated by use of the cranking ap- 

jlq . they obtained for ^^Cr 



proximation 
E{2+) = 0.82 MeV 



15|. In Ref. 



i?4/2 = 2.2, Q(2+) = -34 efm2, and 



B{E2 t) = 1515 e^fm^. The B{E2) value and g(2+) 
are similar to our results, whereas E{2'^) and -E(4^) are 
rather larger. With the DlS-Gogny EDF, there is still a 
clear minimum at the spherical point, and a shallow local 
minimum around /3 = 0.3 16 1. 



In our approach, the time-odd components of the mean 
field are fully included for the calculation of the collec- 
tive masses. Figure EJa) shows the calculated vibra- 
tional masses along the collective coordinate. The col- 
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FIG. 4: (Color online) Collective wave functions of the 
Oi state in Cr isotopes. 



lective masses calculated by use of the cranking approx- 
imation are also shown, where the time-odd components 
are neglected. The cranking masses show a smooth be- 
havior as functions of the deformation /?. This is be- 
cause the cranking mass is evaluated by combination of 
the sum-rule values. In contrast to the cranking masses, 
the LQRPA masses strongly depend on the deformation, 
microscopic structure of the collective mode. And the 
LQRPA masses are larger than the cranking masses. In 
Fig. EJb), the excitation energies squared of the col- 
lective modes are shown. This quantity represents the 
curvature of the collective potential as seen in Eq. (|12|) . 
The two-humped structure is associated with the exis- 
tence of the two minima in the collective potential of 
64Cr. 

Finally, we discuss the shape evolution of the ground 
state. Figure |4] shows the vibrational wave functions of 
the of state. In ^^Cr, the wave function is distributed 
around /? = 0. But the wave functions is not sharply lo- 
calized at the spherical minimum of the potential, which 
reflects the softness of the collective potential against 
the deformation. When two neutrons are added to ^^Cr, 
broadening of the wave function can be seen. When two 
more neutrons are added, the wave function moves to- 
ward a prolately deformed minimum with a large spread- 
ing. From this figure, we can say that ^"^Cr is located 
close to the critical point of the shape-phase transition 
in neuron-rich Cr isotopes and the large-amplitude dy- 
namics plays a dominant role. 

In summary, we have developed a new framework of the 
microscopic model for the large-amplitude collective mo- 
tion based on the nuclear EDF method. The collective 
masses and the potential appearing in the quadrupole 
collective Hamiltonian are evaluated employing the con- 
strained HFB and local QRPA approach, where the time- 
odd components of the mean field are fully taken into ac- 
count. The microscopic Hamiltonian is constructed from 
the Skyrme and the pairing EDFs. 

We applied this new framework to the shape-phase 
transition in neutron-rich Cr isotopes around = 40. 
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The present calculation gives consistent results for the 
low-lying excited states with the observations and the 
other theoretical approaches. Investigating the collective 
wave functions, we reach a conclusion that ^°Cr is located 
close to the critical point of the shape phase transition, 
and the onset of deformation takes place at = 38. The 
large-amplitude dynamics plays a dominant role in the 
shape changes in neutron-rich Cr isotopes. 

As a future work, it would be interesting to compare 
the present framework with the generator-coordinate 
method using the /3-constrained HFB states. And the 
systematic investigation of the properties of low- lying col- 
lective states such as E{2^), -R4/2 and B{E2) in the entire 
region of nuclear chart is quite challenging both in nu- 
clear structure physics and in computational science as 
a large-scale calculation employing the massively-parallel 
supercomputers . 
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